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ABSTRACT 

Context. Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) constitute a special population of 
young neutron stars, which are thought to be magnetars, i.e., neutron stars with super-strong magnetic fields (10 14 — 
10 15 ) G. 

Aims. Assuming that AXPs and SGRs accrete matter from a fallback disk, we attempt to explain the energy-dependent 
pulse profiles exhibited by the AXP 4U 0142+61, as well as its phase-dependent energy spectra. 

Methods. We test the hypothesis that not only the X-ray spectra, but also the energy-dependent pulse profiles of 4U 
0142+61 are produced by accretion along dipole magnetic field lines of strength 10 12 — 10 13 G at the neutron-star 
surface. 

Results. In the fallback disk model, the Thomson optical depth along the accretion funnel is significant and bulk-motion 
Comptonization operates efficiently. This is enhanced by resonant cyclotron scattering. The thus scattered photons 
escape mainly sideways and produce a fan beam, which is detected as a main pulse up to energies of ~ 160 keV. 
The approximately isotropic emission from the stellar surface (soft thermal photons and reflected hard X-ray ones) 
is detected as a secondary pulse. This secondary pulse shows a bump at an energy of about 60 keV, which may be 
interpreted as resonant cyclotron scattering of fan-beam photons at the neutron-star surface. This implies a dipole 
magnetic field strength B « 7 x 10 12 (1 + z) G, where z is the gravitational redshift. 

Conclusions. Our model explains not only the soft and hard X-ray spectra of the AXP 4U 0142+61, but also its energy 
dependent pulse profiles. If our interpretation of the energy dependence of the secondary pulse is correct, the surface 
dipole magnetic field strength is comparable to that of X-ray pulsars. Much like our Sun, the surface multipole magnetic 
field strength may be two orders of magnitude larger, thus allowing for energetic bursts to occur. The accretion process 
is mediated by the dipole field and is not interfering at all with the multipole field. 

Key words, pulsars: individual (4U 0142+61) - X-rays: stars - stars: magnetic fields - accretion disks 



1. Introduction 

Anomalous X-ray pulsars (AXPs) and Soft Gamma Ray 
Repeaters (SGRs) are young neutron stars with X-ray lu- 
minosities much larger than their spin-down power and long 
periods in the range 2 — 12 s. They are widely believed to 
be magnetars deriving their X-ray emission from the decay 
of super-strong magnetic fields (> 10 15 G) (e.g. Duncan & 
Thompson 1992; Thompson & Duncan 1995). At the out- 
set, the magnetar model was developed to explain the gi- 
ant bursts and the large bursts of SGRs, which exceed the 
Eddington limit of neutron-star luminosities by a very large 
factor. If the observed spin-down of these sources is inter- 
preted as the consequence of magnetic dipole braking, the 
resulting polar field strengths are of the order of > 10 15 G, 
in qualitative agreement with what has been inferred from 
the observed luminosities (Kouveliotou et al. 1999). Later 
on, it was discovered that AXPs show short bursts as well, 
though less frequent and less energetic, leading to the gen- 



eral notion that both types of sources are closely related or 
represent even a single class. In the magnetar picture, the 
steady X-ray emission of these sources, which have lumi- 
nosities of typically a few times 10 35 erg s _1 , is thought to 
be caused by a twist of the magnetosphere leading to the 
amplification of the magnetic field and the acceleration of 
particles, which produce the X-ray emission. The twist is 
caused by rotational motions of a crust plate and has a life- 
time of > 1 yr (Beloborodov & Thomson 2007) . For recent 
reviews of the magnetar model see Woods &; Thompson 
(2006) and Mereghetti (2008). 

An alternative energy source for the persistent and tran- 
sient X-ray luminosities of AXPs and SGRs is accretion 
from fallback disks, first proposed, but soon abandoned, 
by van Paradijs et al. (1995), followed by Chatterjee et 
al. (2000) and Alpar (2001). This class of models was de- 
veloped further in a series of papers (Ek§i & Alpar 2003; 
Ertan & Alpar 2003; Ertan & Cheng 2004; Ertan et al. 
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2006; Ertan & Cahskan 2006; Ertan et al. 2007; Ertan & 
Erkut 2008, Ertan et al. 2009), which make the fallback-disk 
idea quite attractive. The fallback-disk model gets support 
from the discovery of IR/optical radiation from two of the 
AXPs, 4U 0142+61 (Wang et al. 2006) and IE 2259+586 
(Kaplan et al. 2009), which has been successfully inter- 
preted as disk emission. Further indirect support comes 
from the discovery of three planets orbiting a neutron star 
(Wolszczan & Frail 1992; Wolszczan 1995). The fallback- 
disk model explains the spin-down of AXPs and SGRs by 
the disk-magnetosphere interaction and requires only "nor- 
mal" neutron-star dipole fields (10 12 — 10 13 G). In addi- 
tion, it is successful in predicting the period clustering of 
AXPs/SGRs in the range of 2 — 12 s. On the other hand, this 
model cannot explain the super-Eddington bursts, which 
are relatively rare. They are attributed to magnetar-type 
activities occurring in local multipole fields (star-spots). 

The main physical questions addressed in the present 
paper are: 

1) What is the source of energy powering the persistent 
emission of these sources? Magnetic-field decay or accre- 
tion? 

2) How large are the polar dipole-field strengths of these 
neutron stars? Are they super-strong (> 10 15 G ) or "nor- 
mal" (10 12 - 10 13 G)? 

The spectra of AXPs and SGRs in the range 0.5 — 10 
keV have been modeled by a superposition of blackbody 
and power-law or, alternatively, double blackbody emission. 
Recently, hard X-ray tails at energies reaching up to more 
than 100 keV have been observed. For a recent review see 
Mereghetti (2008). The luminosities in these tails observed 
at energies > 20 keV are of the same order as the soft X- 
ray luminosities, which is a challenge to any model. Very 
recently, Enoto et al. (2010) published a homogeneous set of 
data on 7 sources showing both soft and hard components, 
as well as their periods and period derivatives. 

Using the bulk-motion and thermal Comptonization 
(BMC/TC) model of Farinelli et al. (2008), Triimper et al. 
(2010) have shown that accretion from a fallback disk can 
produce the hard X-ray (> 2 keV) spectrum of 4U 0142+61. 
The mechanism operates close to the hot thermal mound 
near the bottom of the accretion column, where the free- 
falling matter is stopped. The required seed photons are 
produced in two places (see Fig. 1): a) The photosphere 
of the neutron star, which emits most of the low-energy 
photons in the 0.5 — 2 keV range and b) the hot ther- 
mal mound at the bottom of the accretion column (e.g., 
Basko & Sunyaev 1976; Becker & Wolff 2007), which has a 
higher temperature than the photosphere. In the following, 
we discuss the formation of the different spectral compo- 
nents of 4U 0142+61 and show that not only the observed 
broad-band spectrum, but also the energy- dependent pulse 
profiles (den Hartog et al. 2008) can be explained by the 
BMC/TC model. In particular, we present evidence that 
the magnetic-dipole field at the surface of this neutron star 
is ~ 8 x 10 12 G, consistent with the expectations of the 
fallback-disk model. 



2. The model 

In what follows, we give the basic picture that we have in 
mind, part of which was presented in Triimper et al. (2010). 
In particular, we discuss the geometrical constraints that 
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Fig. 1. Schematic description of the beaming pattern. The 
different wiggly arrows show the mean direction of the re- 
spective component. Blue: The fan beam (main pulse), con- 
sisting of photons with energies between ~ 2 and 160 keV, 
is produced by the BMC/TC process. It is subject to grav- 
itational bending. Green: Fan-beam photons hitting the 
neutron-star photosphere will be reflected (scattered) . Red: 
Soft photons from the photosphere escape, while a small 
fraction contributes as seed photons for the Comptonization 
process. 



the observations of 4U 0142+61 place on the emission re- 
gion and the accretion column and then we discuss the role 
that the dipole magnetic field plays in the formation of the 
X-ray spectrum and the energy- dependent pulse profiles. 

2.1. Geometrical constraints for the emission region and 
optical depths 

Fig. 1 shows a schematic picture of the X-ray emitting re- 
gion. Following Triimper et al. (2010), we assume that the 
hard X-ray component (2 — 160 keV) is produced by ther- 
mal Comptonization (TC, 2 — 10 keV) in the accretion shock 
and the hot mound at the bottom of the accretion funnel 
and bulk-motion Comptonization (BMC, > 10 keV) in the 
accretion flow. 

Despite the fact that the magnetic field plays a signifi- 
cant role in the radiative transfer (see Section 2.2 below), 
it is instructive to look at the Thomson optical depths in- 
volved in this problem. In the free-fall regime, the trans- 
verse Thomson optical depth r t in the direction perpendic- 
ular to the magnetic field, near the neutron-star surface, 
is 

M 

T t = n a T a = nCr T a , (1) 

mpVffiraf, 

where ao is the radius of the accretion funnel at the neutron- 
star radius, no and Vff are the electron density and the free- 
fall velocity at the bottom of the accretion funnel, M is the 
accretion rate, <jt is the Thomson cross section, and m p is 
the proton mass. For M ss 2 x 10 15 g s _1 (corresponding to 
an X-ray luminosity of 3.7 x 10 35 erg s _1 ) and Vff » 0.5 c, 
we get for the transverse optical depth r t ~ 1.7 x 10 4 /ao, 
where ao is measured in cm. In order to have the BMC 
process work efficiently, the BMC region has to be mod- 
erately optically thick in the transverse direction, which 
implies that ao must be of the order of 100 m or less. On 
the other hand, the blackbody radius r^b of the observed 
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soft component (0.5 — 2) keV emitted by the photosphere 
is much larger. It was found to be 8.6 km (White et al. 
1996), 5.4 - 8.3 km (Israel et al. 1999) or 7.2 km (Juett et 
al. 2002). All these radii are normalized to a distance of 3.6 
kpc. 

In the direction parallel to the magnetic held the 
Thomson optical depth of the free-fall zone (the BMC zone) 
is 

rp-^-^l.lxlO 10 /^. (2) 

For ao 100 m, we obtain r p > 100. In this direction the 
Thomson optical depth is huge and the accreting matter 
acts like an impenetrable wall to the outgoing photons. 
Because of this, we can find an upper limit to the trans- 
verse Thomson optical depth r t as follows: 

The hard X-ray spectrum, say from 10 to 100 keV, is 
produced by repeated "head-on" collisions (Mastichiadis & 
Kylahs 1992) between the photons that come out of the 
thermal mound and the infalling electrons, whose Lorentz 
7 is approximately 1.15 near the neutron-star surface. 
The number of such collisions required to upscatter an 
E t ~ 10 keV photon to an E f ~ 100 keV one is N sc ~ 
\og(E f / Ei) / log j 2 ~ 8. On the other hand, the observed 
photon- number spectral index of the hard X-rays is T ss 1 . 
This means that after each head-on collision, roughly 25% 
of the photons escape, while the rest suffer yet one more 
head-on collision with the infalling electrons. 

Let now H be the mean free path of the outgoing ther- 
mally Comptonized photons in the infalling plasma close 
to the neutron-star surface. Then, H « l/(n <JT) — ao/ T t- 
The upscattered photons escape sideways, through a cylin- 
drical surface of height ~ H . The fractional solid angle that 
this surface subtends at the axis of the accretion funnel is 
Sl csc /27r ~ H/a ~ l/r t . Equating this with the escape 
probability of ~ 0.25, we hnd r t ~ 4. 

2.2. The role of the dipole magnetic field 

The Farinelli et al. (2008) BMC/TC mechanism applied 
by Triimper et al. (2010) has been designed for the non- 
magnetic case, since it is using Thomson cross sections. 
Since we assume that the polar magnetic holds of the AXPs 
and SGRs are of the order of 10 12 — 10 13 G, magnetic cross 
sections (e.g. Meszaros 1992) should be used. After the 
pioneering work of Basko & Sunyaev (1976), many peo- 
ple have worked in this area. Becker & Wolff (2005, 2007) 
have advanced this complicated undertaking by an ana- 
lytical approach using angle averaged magnetic cross sec- 
tions and a cylindrical approximation for the accretion col- 
umn. However, their results are not directly applicable to 
AXPs/SGRs, because they have been obtained for sources 
like Her X-l and LMC X-4, which have higher luminosities 
by more than a factor of a hundred. 

The most important effect of the strong magnetic held 
is that the cross section of the extraordinary-mode photons 
is significantly different from the Thomson cross section. At 
photon energies E much smaller than the cyclotron energy 
E c = (h/2ir)(eB/m e c), the cross section becomes much 
smaller than ot for photons propagating parallel to the 
magnetic held, while it becomes very large close to the cy- 
clotron resonance. These effects are softened somewhat by 
the effects of vacuum polarization and spin-flip scattering. 
We expect that the inclusion of the magnetic effects will 



lead to a change of the values of the parameters ao, 7t, and 
t p (Section 2.1). However, the inclusion of these effects will 
not change the main result of Trumper et al. (2010), namely 
the fact that the broad-band (0.5 — 200 keV) spectrum of 4U 
0142+61 can be represented by a single physical model in- 
volving TC and BMC of photons by the infalling electrons. 
Any deviation from the geometry of the TC/BMC model 
of Farinelli et al. (2008) and the inclusion of effects due to 
a strong magnetic field will be absorbed in the values of 
the model parameters and as a result the later should not 
be taken at face value. We expect that the parameters re- 
lated to the seed photon spectrum (temperature and slope 
of the modified blackbody, and the illumination factor) will 
be mostly affected by the change in the geometry discussed 
in Section 2.1 (Fig. 1). 

Extraordinary photons of energy E will be scattered 
quasi-isotropically by electrons when passing sheets in 
which the resonance condition E ss E c — j(h/2ir)(eB/m e c) 
is fulfilled, where 7 is the Lorentz factor of the electrons 
(7 w 1.15 for free-falling electrons near the neutron-star 
surface and 7 = 1 for thermal electrons). For photons 
of energy E, the scattering sheet is located at a radius 
r(E) = R{E c q/ E) 1 ' 3, , where E c q is the cyclotron resonance 
energy at the neutron-star surface. 

The resonant scattering optical depth in the accretion 
column, parallel to the magnetic field, greatly exceeds the 
optical depth for Thomson scattering (Lyutikov & Gavriil 
2006) 

7rc 5 (\ keV\ 

T rcs ~ Z 10 5 — — , 4 

8r e u c \ E c J 

where r e is the classical electron radius, and u> c , E c are 
the cyclotron frequency and cyclotron energy, respectively. 
This means that, for photon energies in the range 0.2 — 160 
keV, cyclotron resonance scattering is the main source of 
opacity. In the "resonance layer", located at radius r(E) w 
R(E c0 /Ey/ 3 , photons of energy E are quasi-isotropically 
redistributed. 

In addition, the scattering by the cyclotron harmonics 
has to be taken into account, which for a given photon 
energy E takes place at radii r(E) « R(jE c0 / 'E) 1 ^ , with 
j = 2,3,4.... 

Since the dipole magnetic field decreases rapidly in the 
vertical direction, the cyclotron scattering dominates over 
a wide range of photon energies and strongly prevents ex- 
traordinary photons from traveling upstream. At the same 
time, these scattering effects enhance the efficiency of the 
bulk-motion Comptonization process. 

For extraordinary mode photons, traveling at small an- 
gles with respect to the magnetic field direction, with en- 
ergies much smaller than the local cyclotron energy, the 
scattering cross section becomes much smaller than the 
Thomson value. This appears to have the opposite effect 
than the one discussed above. However, the extraordinary 
photons traveling upwards will eventually be stopped by cy- 
clotron scattering when passing a resonance layer at higher 
heights. 

On the other hand, the ordinary-mode photons also 
cannot escape in the magnetic field direction, because the 
Thomson optical depth t p is significantly larger than one 
and, in addition, the photons are advected downward by 
the accretion flow. Since t p is significantly larger than r t , 
nearly all the ordinary- mode photons escape sideways, i.e., 
perpendicular to the accretion flow. In view of the above, 
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Fig. 2. Schematic presentation of the energy dependent 
pulse profiles measured by XMM-Newton, RXTE-PCA and 
INTEGRAL-ISGRI. The main pulse (blue phase) ranges 
from the lowest to the highest energies. In our model, this 
component is mainly formed by the fan-beam photons pro- 
duced by the BMC/TC process. The secondary pulse (green 
phase) is produced by the photospheric radiation at low en- 
ergies (< 2 keV) and the fan-beam photons, which hit the 
photosphere and are reflected there. The intensity of this 
pulse falls off quickly between 0.8 and 8 keV. It reappears 
at higher energies and becomes as prominent as the main 
pulse in the 16.3 — 50 keV interval. This figure is a simpli- 
fication of Fig. 7 of den Hartog et al. (2008). 



we conclude that the column will largely emit a fan beam, 
as shown schematically in Fig. 1. 
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Fig. 3. This figure is a simplification of Fig. 8 of den Hartog 
et al. (2008), from which the data are taken. The black 
line represents the INTEGRAL /XMM-Newton total spec- 
trum fit shown in that paper. The spectrum of the main 
pulse is shown in blue. The spectrum of the secondary pulse 
in shown in green. In our model, this component consists 
above ~ 8 keV of fan-beam photons reflected by the pho- 
tosphere. Its spectrum shows a rapid rise up to energies 
of ~ 60 keV and a break at ~ 80 keV. In our model, this 
bump, which can also be seen in the pulse shapes (Fig. 2), is 
interpreted in terms of enhanced reflection due to cyclotron 
resonance scattering. 



As noted by den Hartog et al. (2008), morphology 
changes are ongoing throughout the whole energy range. 
However, the principal structure of the pulse profiles is very 
simple: It mainly consists of two pulses (c.f. Fig. 2). In 
Figs. 2 and 3 we present simplified versions of the energy- 
dependent pulse profiles and the phase-dependent spectra, 
which arc labeled by colors. The main pulse (blue color), 
peaking around phase 0.1, is located at a stable position 
from the lowest to the highest energies and the secondary 
pulse (green color) is located at phases around 0.6 (Fig. 2). 
The phase-dependent spectra (Fig. 3) show the same color. 



2.3. Interpretation of the 4U 0142+61 pulse shapes and 
energy spectra 

A comprehensive and detailed description of the broad- 
band X-ray characteristics of AXP 4U 0142+61 has been 
given by den Hartog et al. (2008). In the following, we 
present a detailed discussion of these results in terms of 
the accretion model. To this end, we use two results of 
the den Hartog et al. (2008) paper, the energy-dependent 
pulse profiles shown in Fig. 2, measured by XMM-Newton, 
RXTE-PCA, and INTEGRAL-ISGRI, and the complemen- 
tary phase-dependent energy spectra shown in Fig. 3. 



2.3.1. Interpretation of the main pulse 

The interpretation of the main pulse is very simple in our 
model. We identify it with the fan beam emitted from the 
region of the hot shock/mound and the infalling matter 
above it (Fig. 1). At the lowest energies (E < 1.7) keV, 
the emission consists of two components, the blackbody-like 
emission from the part of the extended photosphere that 
faces the observer and the low-energy tail of the BMC /TC 
spectrum that is emitted as a fan beam. In the energy inter- 
val 1.7 < E < 10 keV, the spectrum is composed of the tail 
of the photospheric spectrum and the power-law spectrum 
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produced by the thermal Comptonization (TC). At energies 
E > 10 kcV, up to the highest energies beyond 100 keV, the 
emission is due to the bulk-motion Comptonization (BMC). 

2.3.2. Interpretation of the secondary pulse 

In Figs. 2 and 3 we show as green lines the secondary pulse, 
which occurs at phase ~ 0.6. We identify it with the polar 
beam (Fig. 1). As in the case of the main pulse, here also 
we consider three energy bands. 

The intensity of the secondary pulse rises strongly to- 
wards low energies (Fig. 2). An extrapolation of the trend 
shows that at energies below 0.8 keV it becomes stronger 
than the main pulse, which is fully consistent with the 
ASCA data in the 0.5 - 1.7 keV band (Kuiper et al. 2006; 
den Hartog et al. 2008). We interpret this component in 
terms of the blackbody-like spectrum dominating at ener- 
gies below 2 keV and originating in the extended hot pho- 
tosphere surrounding the accretion column, which has a 
rather large extent (r^ ~ 5 — 8 km) . The main difference in 
the angular distribution compared with that of a blackbody 
is that a broad peak occurs in the polar direction, caused 
by extraordinary photons escaping from larger depths (c.f. 
Pavlov et al. 1994). The emission from this area will be com- 
posed of an approximately isotropic component and a nar- 
rower pencil beam of extraordinary-mode photons, which 
escape from the deeper layers along the magnetic field lines. 

In the 1.7 — 10 keV band, the intensity of the secondary 
pulse falls off steeply with energy. In our model, this part 
of the spectrum is produced by the superposition of the 
tail of the blackbody distribution and the power-law spec- 
trum produced by thermal Comptonization (TC). A possi- 
ble third component can be produced by reflection of a frac- 
tion of fan-beamed TC photons, which hits the polar-cap 
area. The reflected photons leave the neutron-star surface 
more or less isotropically, i.e. with an angular distribution 
similar to that of the thermal emission from the photo- 
sphere. Thus, both components are in phase. 

In the 10 — 50 keV band, the secondary pulse becomes 
quite strong and it is actually comparable in strength with 
the main pulse. We note that this pulse extends into the 
50 — 160 keV energy band, where it is visible as a shoulder 
of the main pulse (Fig. 2). But the effect is small due to 
the large width of this energy interval. It can be seen better 
if one looks at the phase-resolved spectra shown in Fig. 3. 
The secondary pulse spectrum shows a broad bump in the 
20 — 80 keV range. Actually, its intensity rises faster than 
that of the main pulse up to an energy of ~ 60 keV, beyond 
which it indicates a drop. We will return to this feature 
in Section 3. In our model, this part of the spectrum is 
produced by the BMC photons that hit the neutron-star 
surface and get reflected (scattered) with some energy loss. 
Since the reflected photons have a more or less isotropic 
distribution in direction, this energy band is in phase with 
the previous two. 

2.3.3. Remarks 

Concluding our interpretation of the energy-dependent 
pulse shapes, we think it is appropriate to give some re- 
marks: 

1. In Fig. 1, we show the swing of the line of sight dur- 
ing half a phase cycle, approximately required by our model 



for 4U 0142+61, which ranges from a direction roughly per- 
pendicular to the accretion column (main pulse) to an off- 
pole direction (secondary pulse). We stress that the phe- 
nomenological appearance of an AXP or SCR in terms of 
energy-dependent pulse shapes or phase-dependent spectra 
will depend crucially on this viewing geometry, which in 
turn depends on the angles between the rotational axis of 
the neutron star and the line of sight on one hand and the 
dipole axis, on the other. 

2. The gravitational bending of the fan beam may play 
an important role for the beaming, depending on the un- 
known mass and radius of the neutron star. 

3. As mentioned above, the hot polar cap has a radius 
r bb ~ 5 - 8 km. The spectral fits of White et al. (1996), 
Israel et al. (1999), and Juet et al. (2002) yield a temper- 
ature of kT w 0.4 keV. However, this is probably a flux- 
averaged temperature, since the temperature will decrease 
from higher values close to the accretion column to lower 
values at larger distances. 

4. The energy-dependent pulse shapes in Fig. 2 repre- 
sent only the pulsed components without the constant com- 
ponents. The pulsed fractions can be found in Fig. 9 of den 
Hartog et al. (2008). They are quite low at soft energies 
(~ 10% at 1 keV), which is qualitatively expected from our 
model, since the photosphere is visible to the observer all 
the time (see Fig. 1). On the other hand, the TC/BMC 
components show larger pulse fractions (~ 25% above 28 
kcV), as anticipated for the dominating fan beam. 

3. An estimate of the magnetic dipole field of AXP 
4U 0142+61? 

We believe that the dipole magnetic field is not only chan- 
neling the accretion flow onto the poles, but also leaves 
its imprint in the observed spectra. As already discussed 
above, the spectrum of the secondary (polar) beam is char- 
acterized by a fast rise (faster than that of the fan beam) 
from ~ 10 keV up to a maximum at ~ 60 keV, followed 
by a decline or a cutoff beyond ~ 70 keV. This behavior 
is expected if the polar beam is produced by scattering 
at the photospheric plasma, whose efficiency rises for en- 
ergies approaching the cyclotron energy. Interpreting the 
spectrum this way, leads to a cyclotron energy of 70 — 80 
keV, corresponding to a polar magnetic field strength of 
B ~ 6 — 7 x 10 12 (1 + z) G, where z is the gravitational 
redshift. Clearly, observations with higher photon statistics 
would be necessary to confirm this conclusion. 

4. Heating of the hot polar cap 

The photosphere responsible for the soft component has a 
flux-averaged temperature of ~ 0.4 kcV and a blackbody 
radius r b b ~ 5 — 8 km. In our model, there are three possible 
ways to heat the polar cap: 

1. By absorption of a fraction of the fan-beam photons 
hitting the neutron-star surface. This process will be en- 
hanced by the gravitational bending of the fan beam to- 
wards the neutron-star surface, which depends on the grav- 
ity of the neutron star, i.e., its mass over its radius. If this 
were the main heat source of the polar cap, the power of 
the re-emitted thermal component and the reflected (scat- 
tered) component would be equal to the power of the fan 
beam impinging on the polar cap. 
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2. A second heat source may be that part of the accre- 
tional power which is conducted into the stellar interior and 
redistributed over a larger area, depending on the magnetic- 
field configuration in the interior and the crust. 

3. Another contribution to the heating of the polar cap 
could be provided by conduction of heat from the hot inte- 
rior of the relatively young neutron star. 

In view of our limited quantitative understanding of all 
properties and processes, which play a role in this context, 
it would be premature to favor a specific combination of 
these processes. 

5. Summary and conclusions 

We have shown in a semi-quantitative way that bulk- 
motion/thermal Comptonization in an accretion column, 
formed by a dipole magnetic field of strength ~ 10 13 G, 
describes well not only the soft and hard X-ray spectra 
but also the phase-dependent energy spectra of AXP 4U 
0142+61. The energy- dependent pulse profiles and their 
constancy over long periods of time constitute significant 
observational constraints for the proposed models. Our 
model explains them in a simple and natural way by the 
formation of two "beams" , one perpendicular to the ac- 
cretion column (main pulse) and one roughly parallel to 
it (secondary pulse). AXP 4U 0142+61 is one of seven 
sources presently known to have hard X-ray tails (Enoto 
et al. 2010). From the similarity of their broad-band spec- 
tra, it is no surprise that accretion with thermal and bulk- 
motion Comptonization works for all of them (Zezas et al., 
in preparation). 

On the other hand, the magnetar model of a twisted 
magnetic field (Thompson et al. 2000), where the hard X- 
ray emission originates far up in the magnctosphere, seems 
to have problems in explaining the energy-dependent pulse 
profiles (Fernandez & Thompson 2007) Another problem 
for the magnetic-twist model is the long term stability of 
the pulse shape and phases. This would require that the 
rotation of the same crustal plate, located at a fixed po- 
sition on the neutron star, is responsible for the quiescent 
emission, which looks quite artificial. Furthermore, the P 
of SGR 0418+5729 was determined (Rea et al. 2010) to be 
less than 6 x 10~ 15 s s _1 , which in the magnetar model 
implies a magnetic dipole field of < 7.5 x 10 12 G! 

In this paper we propose that the spectral bump of the 
secondary pulse (which corresponds to the polar beam) is 
caused by cyclotron-resonance reflection of the fan beam 
by the magnetized photosphere. This leads to a magnetic- 
field strength of 7 x 10 12 (1 + z) G, which falls in the range 
of dipole magnetic fields (10 12 — 10 13 G) attributed to 
AXPs/SGRs in the framework of the fallback-disk scenario. 

On the other hand, we stress that the few giant bursts 
observed from SGRs cannot be explained by any accre- 
tion process. The same is true for bursts with large super- 
Eddington luminosities. They are most likely caused by 
crustal shifts of plates carrying super-strong magnetic fields 
(>10 15 G), as discussed in the classical magnetar litera- 
ture (e.g. Thomson & Duncan 1995). Our analysis suggests 
that these events do not take place in the dipole field, but 
in localized multipole fields. This situation is qualitatively 
similar to that of the Sun, which shows flare activities in 
sunspot fields, which are larger than the solar dipole field 
by at least two orders of magnitude. 



In conclusion, our model explains naturally the X- 
ray spectra and the energy-dependent pulse profiles of 4U 
0142+61. Furthermore, it allows this source or any other 
AXP or SGR to emit big or giant bursts if the neutron 
star involved sustains multipole fields with strengths in the 
10 14 - 10 15 G range. 
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